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In this article we propose the calculation of the unconditional Wiener measure functional integral 
00 . with a term of the fourth order in the exponent by an alternative method as in the conventional 

' perturbative approach. In contrast to the conventional perturbation theory, we expand into power 

I series the term linear in the integration variable in the exponent. In such a case we can profit from 

the representation of the integral in question by the parabolic cylinder functions. We show that in 
' such a case the series expansions are uniformly convergent and we find recurrence relations for the 

O I Wiener functional integral in the A'' - dimensional approximation. In continuum limit we find that 

, the generalized Gelfand - Yaglom differential equation with solution yields the desired functional 

integral (similarly as the standard Gelfand - Yaglom differential equation yields the functional 
integral for linear harmonic oscillator). 



oo 



I , I. INTRODUCTION 

^ • We will define the continuum functional integral as the limit of a finite dimensional integral. This finite dimensional 
integral is derived from the continuum one by the time-slicing method. We avoid so problems with continuum integral 
measure because we consider continuum limit of the result of a finite dimensional integral. In order to evaluate finite 
dimensional integrals we must first solve the problem of the calculation of one dimensional integral: 



+ 00 

Ii= J dx exp{-{Ax^ + Bx^ + Cx)} , (1) 



in 
cn 

where Re ^ > 0. The exact analytical result for this integral is not known yet, therefore one uses approximative 
methods of calculation. 

OO , The usual perturbative approach is based on Taylor's decomposition of the fourth order term with consecutive re- 
■ placements of the order of integration and summation: 

oo 

X ■ ^^4^ / dx x^" exp{-(B.x2 + Cx)} (2) 



n=0 



The above integrals can be calculated, but their sum is divergent. 

However, Ii = Ii{A, C) is an entire function for any complex values of B and C, since there exist all integrals 

ff^d'ShiA, B, C) ^ (-1)"+'" J dx exp{-iAx^ + Bx^ + Cx)} , ReA > 0. 

— oo 

Consequently, the power series expansions of Ii — Ii{A, B. C) in C and/or B have an infinite radius of convergence 
(and in particular they are uniformly convergent on any compact set of values of C and/or B). We shall frequently 
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use the power series expansion in C: 



h = Y / dx exp{-(Aa;4 + Bx^)} 

n=0 J!^ 



(3) 



The similar idea of the first order term expansions in the fourth order action was used by Tuszyhski et al. [l| for an 
evaluations of a non-Gaussian models for critical fluctuations of the Landau - Ginsburg model of phase transitions. 

The integral in ([3]) can be expressed in terms of the parabolic cylinder function D^{z), (see, for instance, [^). For 
n odd, due to symmetry of the integrand the integrals are zero, for n even, n = 2m we have: 



V4 



(V2l)"'+i/ 



— r(m + l/2) D_,„_i/2(z) 



dy 2/' 



1-1/2 



exp{-Ay^ - By} 



where 



Explicitly, for the Eq.(l3]) we have: 



r(i/2) ^ (0 



1/4 



^ ' m=0 



B 



-D_ 



1/2(2) , C = 



C2 



4V2A 



(4) 



This sum is convergent for any values of C, B and A positive. 

The convergence of the infinite series in Eq. ^ can be shown as follows. For \z\ finite, \z\ < and | arg{—i') |< 
7r/2 and if | |— > 00, the following asymptotic relation is valid 



D^z) = -i= exp ^(In - 1) - z 

The m term of the sum in Eq. possesses the asymptotic: 

1 



1 + 



■ exp 



(m+^1/2) ^^^ + 1/2) - 1) - v/(m + 1/2) z + mln^ 



(5) 



(6) 



This means, applying the Bolzano-Cauchy criterium that the sum in Eq. Q is not only absolutely, but uniformly 
convergent for the finite values of the constants of the integral (H]) . 

II. EVALUATION OF THE FUNCTIONAL INTEGRAL BY TIME SLICING METHOD 



We suppose that Gaussian integration over momenta is done. Our aim is to evaluate the continuum unconditional 
Wiener measure functional integral: 



where the continuum action contains the fourth-order term: 







(7) 



The functional integral Z is defined by a limiting procedure from the finite dimensional integral Zjy, obtained from 



the continuum integral, when the infinite measure [I?(y9(a;)] is replaced by the finite dimensional measure Y\^^i dLpi{. 
i: 



-■N 




exp 



N 

E 

1=1 



A 



c/2 



- 
A 



4 



(8) 
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where A = /3/iV. The unconditional measure integration is characterized by integration over variable ip^. This is the 
only difference from conditional measure integration, when the ip^ variable is fixed. The detailed discussion of the 
conditional Wiener measure case is technically more involved and was not in the program of this article. We evaluated 
this case in a rather simplified form in the Appendix 4 of our article I . We included a brief discussion of this case 
in conclusions. The continuum unconditional Wiener measure functional integral is defined by the formal limit: 



lim Zjv 

JV^oo 



To evaluate this limit we follow the idea of the Gelfand-Yaglom proof of the functional integral for the harmonic 
oscillator Q, based on the iterative procedure for the finite dimensional representation of the functional integral. The 
idea of N dimensional integration of ([S]) is explained in Appendix A, we quote here the result: 



Zat = [27r(l + 6AVc)] ' [2n{l/2 + bA^ /c)] 



-1/2 



N 



E n 

^1 i^N — 1—0 4 — 1 



2ki 



-r(fc,_i + h + 1/2)1?- 



fei-i-fci-i 



/2 (Zt) 



(9) 

where ko = = 0, = ^2 = ■ ■ ■ = = ^ = (1 + bAyc)-\ ^^-i = (1 + 6AVc)-i/2(l/2 + bAyc)-^/^ = , 
= Z2 = ■ • • = ZN-1 = z = c{l + 6AVc)/v/2aA3, zn = c{l/2 + bA'^ /c)/ ^/2^. 



III. GENERALIZED GELFAND ~ YAGLOM EQUATION (GGYE) 

Let us rewrite the result for the N dimensional integral ([9]) in the form: 

" N 

Zn= Y[2{l + bA^/c)LJ, S 



N 



with 



00 N 

E n 

^1 1 ■ ■ "T^iv — 1— i—1 



(6) 



2ki 



r-i(i/2)v^ r(fc,_i + h + i/2)v_k,_,-k,-i/2 {z^) 



(10) 



(11) 



where the constants and symbols in the above relation are connected to the constants of the model by the relations: 
LOi = 1— A^/wi-i, Wo = (l/2-|-6A^/c)/(l+6A^/c), A = 2(i+fcA^/c) • We prove the above form (fTU|) of the N dimensional 
integral © later, now we use Eq. (|10p for explanation of the Gelfand - Yaglom procedure of the construction of the 
difference equation. This difference equation is converted to differential equation in the continuum limit N 00 . 
Let us define functions by: 



n 2(l + 6AVc)c^, 

i=0 

''k 



^fe = — F-2 (12) 



The function Fk is defined from the relation of the N dimensional integral pi7|) with A fixed and the variable N 
The quantity given in Eq. ^ is related to F/v as follows: 




The aim of the Gelfand-Yaglom construction is to find the continuum limit of the difference equation for the function 
Fk- Solution of this differential equation is connected to the continuum functional integral by: 



where /3 is the upper bound of the time interval in the action ([7|. 
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The idea of the GGYE construction is based on the recurrence form for the factor uji . We replace LUi by the functions 
F^.±i and Sk,k±i ■ Pedagogical descriptions of this procedure can be found in Appendix B, there ia a proof of the 
lemma: 

Lemma. Let Fk be the function defined by: 



n 2(l + 6AVc)a 



Fk 



i=0 



Si 



(13) 



with LOi defined by recurrence relation: 



and ojQ = (1/2 + 6A^/c)/(l + 6A^/c), A — 2{i+bA^/c) ■ "^^^ constants b, c, A are parameters of the model. 
Let in continuum limit the following condition is valid: 



lim (A Oi + A^ O2) = 



(14) 



where 

c 



Sk+i — Sk 
AS, 



fc+i 



, {Fk — Fk-i) 
' A 



Sk+i — Sk 



Sk+i — Sk 
ASk+1 



-Fk-i 



Sk+i — '^Sk + Sk- 
A^Sk+1 



ASk+i 
2 



-4Ffe_i 



^fe+i — Sk\ f Sk+i — 2Sk + Sk-i 
ASk + 1 J \ A^Sk+i 



(15) 



Let in the limit A — )■ 0, the functions F(t) and S{t) are the continuum limit of the function Fk and Sk, when we 
define: lim^^g k.A = t . 

Then F{t) is the solution of the differential equation: 



_f(,)+4(-F(r) 



d 



Ot 



2b 



d 



- In5(r) = F(r) - - 2— ln5(r) ~ 4(- ln5(r)) 



(9r2 



with initial conditions 



m 



1 



OT 



52(0) 

d 



e(0,/3) (16) 



(17) 



dr V52(t) 



T=0 



The nontrivial dynamics is hidden in the function S{t). 

A'^oie/Whcn S{t) is known exactly the above equation can be simplified by the substitution: 



F{t) 



S^r) 



For the new variable 2/(t) we find a simple equation: 



^y(r)=2/(r)(- 



(18) 



accompanied by initial conditions: 



y(0) = l ''^^^ 



dr 



= 0. 

T = 



Thus, in the case when function S{t) is known exactly, the problem of the functional integral calculation is trivial. 

Problems arise in situations, when S{t) is known approximately, as a result of a perturbative approach. Below we 
define a reasonable approximation of S{t) valid in the proximity of t = 0. But for finite (large) t — [3 this asymptotic 
expansion of S{t) does not have to be valid. However, in our case, the development of the function F{t) from r = to 
T = /? is controlled by a differential equation. Approximative knowledge of the function S{t) leads to a more reliable 
result for F(/?) as the solution of Eq. This philosophy of the calculation corresponds to ideas of evaluation of 

physical quantities be the renormalization group approach. 
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IV. EVALUATION OF THE FUNCTION 



The exact result of N— dimensional integration ([9]) contains summations of products of parabolic cylinder functions. 
We were not able to find a simple formula fur such a summation (although, parabolic cylinder functions belong to 
the representation of the group of the upper-triangle matrices, therefore, due to the group theoretical background 
we believe to the simplification of the product of two such functions) . In what follows, we explain our approach to 
provide the summation over indexes ki. For a single index we are dealing with the sum of the series: 

' "fc. = J2kY\^^^'-^ + l/2)I?-fc._i-fe.-i/2 + h+1 + l/2)I?_fc,_fc^^,_i/2 (z) (19) 



The series (|T9|) is uniformly convergent. For asymptotic values of ki following the Stirling formula for logarithm of 
gamma functions and the asymptotic relation ([S]) for parabolic cylinder functions we find 

Inafe^ = -(fc, - ^ ~ '^ )lnfc, + 2h{\n{^z) - In 2 + 1/2) - z + (fc,_i + ki+i + l)lnz 

+ l/21n7r + zV2 + o(fc7^''^) (20) 
The leading term of the above relation is 

ln(afcj ^ -fciln(fci). 



and 



(^z)2fe» (%/^z)'=-i+'='+i 
ki\ exp {2^/ki z) 



This asymptotic behavior of afc. is sufficient for a proof of the uniform convergence of the series. The convergence 
criteria are fulfilled for an arbitrary group of sums over k^ index in the result ^ for N— dimensional integral. 

In fact in ^ we have an iV-tuple sum of two index quantities (closely related to a^i ) ■ Let us for simplicity discuss 
the N tuple sum of the two index quantities: 

CX) 

^ ^ CKfeo,fclQ^fcl,fc2'^fc2,fc3 ■ ■ ■ '-'^fcjV-l,fcNQ^fcjV,fc]V + lI 

fci, AT —0 

afc._i,fc. = ^g^r(Vi+fc^ + l/2)I?-fc,_,_fc,_i/2 (2.) (21) 

If each individual sum over ki exists and is finite, we divide the sum over the index ki to a finite principal sum and 
a remainder ei(A'o), which can be made as small as possible by suitable selection of the upper summation limit Kq^ 
comparing it to the principal part of the sum: 

oo Kq oo K'o 

afc,_i,fcia/ci,fci+i = ^ afc,_i,feiafc,,fei+i + ^ Q!fci_i,fc.afc,,fc.+i = {'^ + £i{Ko)) ^ aki_i,kiaki,ki+i (22) 

ki=0 ki=0 ki=Ko + l ki=0 

Let us define the " principal sum" E(m, N — m) of (j2ip as follows: the first m summations run to infinity and the 
last N — m summations run to Kq : 

oo 

Y.{m,N-m)^ ^ Ofci.fcjafcj.fca ' • ' a/c„_i,fc„&fc„ , (23) 

ki,---,k,„=a 

where 

bk^= X! ak„,^k^^tak,r,^i,k^+2- ■ ■ <^kM,kN+i (24) 

A:^+i,---,A:jv— 
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In this notation, the sum ((2T|) is S(iV, 0), our aim is to estimate this sum by its principal sum S(0;iV), when all 
summations run over finite range. Performing in I](m, N — m) the sum over fc^, we obtain: 

i;(m, TV - m) = (1 + e„OS(rn - 1, TV - m + 1) . 

We have the inequality: 

S](0, N) < S(7V, 0) < (1 + £)^S(0, N) . 
For any fixed N we choose ii'o so that e = max£i(i^o)- 

Following this discussion, the finite dimensional integral will converge if the upper bound e to remainders approaches 
to zero as 

iV-i-«, 6* > 0. 

This will guarantee the convergence of S(0, N) to the continuum integral. 

In what follows, we describe the idea of the evaluation of "principal sum" for Eq. ^ and how to estimate the 
remainder to this leading term. The procedure of the summation consists of the following steps. At first, we represent 
one of the parabolic cylinder functions 'D^m^i/2{z) by Poincare - type expansion Q, valid for real index and positive 
argument of the function, under the assumption that the index of the function is finite and the argument is going to 
infinity. For the dimension of the integral sufficiently great this is consistent, because z ^ N^^"^. We have: 



where ej{m, z) is the remainder of the Poincare - type expansion of the V function. For Poincare - type expansion 
the upper bound of remainder was calculated by Olver Q . We use the improved upper bound evaluated by Temme 
0. The upper bound for remainder in definition (|25p reads: 



where 

2m 



16 



y2 



2z^J (z2-2m)2 



The estimate p6p is valid for 2y/rn < z [8]. Before insertion of the Poincare - type expansion (|25p into (fT9| for one 
of the V, we divide the sum over ki into two parts. One, over finite ki, where Poincare - type expansion is correct 
and the second, the remainder, small compared to the first part due to uniform convergence: 

E -1(2 E 7|i^^(^'-i + ''^ + ^/2)r(/c, + fc.+i + 2j + l/2)I?_,._,_,._i/2 (z) + (27) 
+ E 7|^r(/c,_i + fc, + l/2)r(/c, + fc,+i + l/2)I?_fe,_,_fc,_i/2 (z) ejih + h+i,z) + 

ki=Q ^ 

CO /y-\2ki 

+ E tIfW^^^'-I +h + l/2)r(fc, + h+i + l/2)V_k^_^_k,-l/2 iz)V_k,-k^^,-l/2 (z), 

where 2\/Kq < z. In the next steep we extend the summation in the first, ("leading") term up to infinity by adding 
and subtracting the terms, allowing to sum over the index ki according to the relation[3j: 



^.V4^ ^i'^^-.-fc(^) = e(^-*)^/4i?.,(x-0 . (28) 

fc=0 

We summed up the product of the two functions V , the result is the function V with a new argument. The pedagogical 
description of this procedure can be found in article I Q . We show the idea of the evaluation of the leading term in 
Appendix C. 



7 



We have shown, that the leading term for the N dimensional integral is of the form: 



where the symbol {N)^^ is defined by the recurrence relation 



^ (AQXl^y-^ ^ (^Qa-.Qa-i)^^- (A - 1)11, ( 5; ) iQUr-' ,pe<0,2,> 



(30) 



where the recurrence in {^)\^^_p starts from: 



^0 



Putting A = A^, and p = 2/i in Eg. ([50]) we obtain the relation for the leading part of the N— dimensional integral, 
which in the continuum limit coincides with the functional integral searched, if the remainder of the leading part 
in the continuum limit disappears. The sum over index /i in Eq. (j29p is an asymptotic sum, therefore the upper 
summation limit must be taken symbolically. The estimate of the remainder to Zf^li is discussed in Appendix D. 
We have shown, that this remainder can be made smaller than N^^^^ , 9 > 0, which is the necessary condition for 
the vanishing remainder in the continuum limit N ^ 00. 

V. SOLUTION OF THE RECURRENCE RELATION FOR 



(AQaQa-i)"''^ (A)^'' 



To evaluate the TV dimensional integral, we solve the recurrence relation (|30p for arbitrary value fi — d: 

)2d _ 
2d~p ~ 



The right hand side of the equation is {2d, p)-th matrix element of the product of three matrices. For fixed d on the 
left hand side of the equation, we read only the d-th column of a matrix, which is recurrently tied to the matrix in the 
center of the product on the left hand side. We define an auxiliary matrix X'^(A) by the following matrix equation: 



^L(A) = E E A^,/A-l)C,^,,(A-l)Ml^(A-l). 

3=0 A=[i±i] 

To evaluate the matrix ^(A), we must calculate X''(A) for all dimensions d up to fi, for each dimension to extract 
d-th column of matrix X'^(A) and to compose from these columns the matrix C^(A). We define such linear operation 
as follows: 

1. Let A'' and M"^' are the matrices of the dimensions (2/i + l)(2/i + 1) and (/i + + 1) respectively, is the 
matrix of dimensions (2/i+ + 1). These matrices possess nonzero main minors of the dimensions {2d+ l){2d+ 1), 
(d + l){d + 1), and {2d + l){d + 1), respectively. The definition of the matrices A''(A) and M''{A) is in Appendix E. 

2. Matrix X'' is the one column matrix defined by the relation: 

X'*(A) = X''(A)*P'*, 

where P'' is the projector of the d-th column of the matrix X''(A) into d-th column of the matrix X'^(A). P'' is a 
matrix with a single nonzero term 
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3. The matrix X'^(A) is defined by relation: 

X'^(A) = A'^(A - 1) * C^(A - 1) * M^(A - 1). 

4. Then, for C''(A) we have the result: 

d 

C^iA) = A"*-'^ (A - 1) * C^-"^ (A - 1) * M'*-*'^ (A - 1). 

IA=0 

5. After evaluation of the full recurrence we find: 

d d-iA rf-iA — iA-i ia 

*A=OjA-l=0 12=0 

l^d-iA (A - 1) ^d-iA-.A-l (A - 2) * • • • * . 



d-JA-jA-l 12 



{!)}*€ 



IK—lK-l 12 I 



1) * ■ ■ ■ * 



-i(A-2)*M'^-*'^(A-l)} 



(31) 



d — lA — iA- 1 ^2 



(1)* 



The evaluations of multiple products of the matrices is given in Appendix E. Remember that for the function iSa 
defined in Eq. P5|) the only important matrix element is |C(A)^^}2^ ,^ then we find the result [9|: 



/i /J — lA M~*A— iA-1 is 

E 

12=0 



-E E 

JA=0 JA-1=0 



{C(A)2^} 

2/2 /2 /.T ^ .M ^i^i„J L-i . \ ( A 



i=0 A=|I+i 



^m=2 



/a 

/a-] 



4(/3-/2) ^4(/a-/a-i) 

• ■ ■ ^;a-i 



1 



2/i-2/„ 



+ 



6 • ■ • Ca-] 



A-l 



(a;iC„^l) 



-^2 \ „2A ^4A^4(/2-A) 



A 



^2A-7 Qo 



We introduced the abbreviation 



I J = /i - (ij + ^ Vih). 



are independent variables and is the differential operator given as: 

= 3/4a| + iidl + 

The asymptotic decomposition of the function 5a reads: 

Ai=0 ^ ^ ' 



(32) 



(33) 



since {C(A)2^}2^ = (A)o'^- The quantity z^A-^ is finite in the continumn limit A ^ 00 and the factor A^^ ensures 
that only the leading term of {C(A)2^}^^ survives the continuum limit. Due to the analytic form for {C(A)2'^}^^ 
we can express 5a in the continuum A — > limit as well as in asymptotic ^ 00 limit. The analytical evaluations 
of the lower /x terms and the relation for the asymptotic terms is done in details in our article P . We evaluated the 
analytical result for the first three terms (i.e. /i = 1, 2, 3). As an illustration we give the first of them: 



{C(A)2} 



2,2 



(2^)1 2^ ^ 



j(o,0;2,l;&fc/6fe_i) 



k=2 
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1.0 r 
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0.6 



7T 
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-20 
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20 



Fig. 1: b dependence of the continuum function S{a,b,c,T) for fixed values a — 0.1, c — 0.5, r — 1. The first three nontrivial 
terms of the asymptotic series (|33p were used. 



where 



MIN 



2ij - I 



and 



2tj~p J {n-2i,~p)\ \h,-i^ 



Qj+iQj 

The continuum hmit corresponds to the prescription: 



k.A ^ X, A. A — > T, 



A . 



k=2 



When A ^ A^, then t ^ (3, where [3 is the the constant of the model, and N ~ 



In the continuum hmit we obtain; 



)fe ^ 2 cosh(7a;) 



bk ^ -J — (tanh(7r) — tanh(7a:)) 
A7 



where 7 — y^2b/c, b and c are parameters of the model. The continuum limit of the relation 
S{a, b, c, t). We show the first nontrivial term (/i = 1) of the three evaluated now: 



will be called 



{C^(a, b, c, t)}^ 2 = [37Ttanh^(7r) + tanh(7r) - 77 



(34) 



For the calculated higher terms we have analytical formulas also as results of algebraic evaluation by 
Mathematica[lo| . The continuum function S{a,b, c,t) for the first three nontrivial contributions is shown in Fig. 
1. 

The corresponding term for the Gelfand-Yaglom equation, —2d^\n{S {a, b, c, t)) — 4(9^ ln(5(a, 6, c, r)))^, is shown 
in Fig. 2. 
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2 3 
tau 



Fig. 2: r dependence of the continuum function —2d^ \n{S{a, h, c, r)) — 4(9^ ln(S'(a, b, c, r)))^ for fixed a = 0.l,b — 5,c = 0.5. 
The first three nontrivial terms of the asymptotic series (|33p were used. 



In the limit /i — )■ cx3, A fixed the terms p2p are divergent. For the leading divergent term for fi oo of the 
asymptotic series for 5a ((33|) we find [9|: 



2/j 



(35) 



The series of this form is an asymptotic expansion of the parabolic cylinder function of the index —1/2 and the 
argument 



/ tanh (7r) 



V C7 

For small /i the terms of the decomposition don't fulfil such decomposition. 



VI. CONCLUSIONS 



We presented an analytical method of evaluation of the unconditional Wiener measure functional integral with a 
fourth order term in the action. No simple analytical form of such an integral is known, perturbative methods of 
evaluation are needed. Instead of a standard perturbative procedure we expand the linear term of the action. Such 
integral is an entire function of all remaining parameters with infinite radii of convergence. We find an analytical 
result for the functional integral in the form of the solution of the "generalized Gelfand-Yaglom" (GGY) equation for 
the case of an anharmonic oscillator with positive coupling and no positivity requirement for the quadratic term. We 
calculated the asymptotic solution of the GGY equation up to third order in the coupling constant. 

The same method of evaluation could be applied for the case of the functional integral with the conditional Wiener 
measure for anharmonic oscillator, with more interesting physical results. The approach based on the GGY equation 
can be applied too. However in this case the method is technically rather involved. An evaluation for a simplified case 
when the endpoints are fixed to zero is presented in our article I [5| ■ In this case the functional integral is a propagator 
with coinciding endpoints well known as Moeler's formula Q in the case of harmonic oscillator. We evaluated the 
correction to Moeler's formula S{P) up to the first nontrivial term and we fomid: 



Si(3) = 1 - 



3a 



320273 



-3coth(7/3) + 27/3 



coth^(7/3) 



1 



2sinh^(7/3) 



(36) 



A more complete evaluation of the functional integral with conditional Wiener measure and a more systematic study 
of the anharmonic oscillator is in progress. 

When the frequency 6(r) and the coupling constant air) are time dependent we obtain for the anharmonic oscillator 
by the time slicing method described in the article an equation which in the continuum limit reads: 
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26(7 



2^ In Sir) 



4(-ln5(r) 



(37) 



where S{t) can be evaluated by a method similar to that of Section IV. The time dependence of the functions 
a(r), 6(r), instead of constants a, 5, does not complicate summations over the index in individual time-slice intervals. 
The nonlocal character of the result for TV— dimensional integral, represented by the dependence of parabolic cylinder 
functions on two summation indexes is represented by the function i^^, introduced in the evaluation procedure: 

(1 + bN-,n-lAyc){l + bN-,nAyc)£_%_„^_^ = 1 

After the full recurrence procedure of the evaluation of the N— dimensional integral, we obtain the result: 



1 



TCUt _ J I I 



s 



(38) 



where Sn-i is a function with structure unimportant for this moment. The function tOi is defined by the recurrence 
relation: 



UJi = 1 



We obtained Eq. ([57]) without detailed evaluation of the function S{t), by the same method as in Section ///. By 
substitution 



we convert ([37l) to selfadjoint [ll| equation: 



F{r) 



y"{r) = 



^ v{t) 
26(t 



2/(t)- 



(39) 



In the case of a linear harmonic oscillator with time dependent frequency, we obtain for the function F(t), leading 
to the inverse square root of the functional integral, the simple second order selfadjoint equation discussed for the 
first time by Lewis |12l|. The time evolution of QM systems of a broader class of time dependent Hamiltonians was 
discussed by Samaj [iJl- Equation ([57]) belongs to the class of general linear second order differential equations 
discussed exhaustively by Kamke [ll|. Equation pQ]) corresponds to its reduced normal form with invariant 

. 26(r) 



In Kamke [ll|, the variable substitution and function replacement leading to the differential equation of the second 
order with constant invariant are proven. Therefore, as in the case of the harmonic oscillator, in the case of anharmonic 
oscillator there exists a possibility to convert the problem with time dependent frequency and coupling constant to 
the constant coefficient linear equation of the second order problem. 

Functional integral methods play an important role in quantum mechanics. One has the tools to evaluate the 
mean values of observables without a necessity to solve equations of motion. For example, the conditional measure 
functional integral in quantum mechanics represents particle propagation. Its integrated form over final positions 
describes unconditional functional integral in quantum mechanics. In statistical physics, the same quantity, after 
Wick rotation it — s- — r, represents the partition function. In fact, all such functional integrals offer important 
information about physical quantities, e.g. spectrum of Hamiltonian, mean values of observables, etc. 

The anharmonic oscillator can be considered as 0^ theory in (1 + 0) dimensions and can be regarded as a toy 
model for understanding of QCD, as the anharmonic oscillator [14j studied in the perturbative approach was. In our 
description we have a possibility to study the case with positive mass, corresponding to the anharmonic oscillator 
and, what is more interesting, the case with negative frequency. It will be interesting to analyze the analog of the 
GGY equation in field theory also. 
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APPENDIX A: EVALUATION OF THE iV DIMENSIONAL INTEGRAL 



A pedagogical evaluation of the integral is done in our article I [5], here we recall the most important steps. We 
are going to evaluate the N— dimensional integral defined by the relation: 



+00 



-■N 



N 

n 

i=l 



N 



2-kA 



exp ■ 



i=l 



c/2 



- ■^g^-l 
A 



(Al) 



First, we rewrite the sum in exponential function in a form convenient for consecutive integrations: 

foA2 , c 



A 

+ ••• 

c feA^ c 
+ ^o^ipj + — (1 + — _ —^^^p.^^ 



(A2) 



+ Aa^% + 7-1/2+ )ip% . 

A c 



We expand the exponential factor containing terms linear in the integration variable into Taylor's series. Using the 
integration formula 0: 



^ exp{-px^ - qx) dx = r(a)(2p)-°/2 exp (^^^ 



we find for integration over the variable ipi 

1 



E 



A(l+6AVc) 



(^2)''^ r(fci + i/2)i?_fc^_i/2W , 



v/2^(l + 6AVc) (2/ci)! 
where we used the notation: 



c(l + bAVc) 
^2^ ■ 

The term (1^92)^'^^ in Eq. (jA4p will play active role in the integration over the variable ip2- 



^2 — 



(^2)^''' exp<^-Aa(^^ - — (1 + )<y5^ + -r</'2¥'3 

27rA I A C A 



(A3) 



(A4) 



(A5) 



(A6) 



Taking both integration steps together we have: 

2 ^ c2ki 



Z-2 - 



1 



vV24l + 6AVc)j 
where we used a new symbol: 



A(l+6AVc) 



(2fci)!(2A:2; 



r(A:i + l/2)I?_fe,_i/2(z) r(A:i + fcz + l/2)P_fc,_fe,_i/2(z) , 

(A7) 

1 



(l + 6A2/c) 



Eq. (jA7p is the result of both ip\ and integrations. We can repeat this procedure for other integration variables 
• • • , <y9jv-i. For the integration over variable (^^r one has no linear term in the exponent, therefore we don't expand 
anything and this last integration will not add summation over the index kN to the final formula. We have: 



+00 



C2'l0C 



(^^f^-- exp|-Aa^4^ - + 

27rA I A C 



(A8) 
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and the result is: 



710C 
-■N 



(c(l/2 + &AVc)/A)'="^^ 



^27r(l/2 + feAVc) 
Remember the difference in definitions of zjsi and z.^, i = 1, 2, • • • , iV — 1: 

c(l/2 + 6AVc) 



T{kM-i + l/2)P_fc„_^_i/2(zw) . (A9) 



ZN 



V2^ 



and also the term 



which modify the definition of: 



(c(1/2 + 6AVc)/A)'"-^ 



1 1 



IN-l 



(l + 5AVc) Y (l/2 + 6A7c) 
For the N— dimensional integral we obtain finally the exact result: 



1/9 ^ ^ 

2:Ar= [2^(1 + 6AVc)] ~ [27r(l/2 + 6AVc)] ' ^ J] 



^^^^ ■r(fc,_i + fc,; + l/2)2?_fc,_,_fc,_i/2 (2,) 



where ko = = 0, and = ^2 = ■ • ■ = = also ^tv — 1 , and = 22 = • ■ ■ = 2:a^-i = ^ 



(AlO) 



APPENDIX B: PROOF OF THE GENERALIZED GELFAND-YAGLOM DIFFERENTIAL EQUATION 

We rewrite the N— dimensional integral given in Eqs. (fTO|) . (jlip as follows: 

5jv-i(A) 



Zn — 



N-l 

I n 2(l + 6A2/c)w, 

1=0 



where 



1/2 + 6AVC ^ 1 



= 1 , LUo = , , ■ , , — , A 



Ui-i' " 1 + &A2/C ' 2(l + 6A2/c)' 
The value of the functional integral in the continuum limit is formally defined by 



Z = lim Zjv 

N^oo 



Let us define the function 

k 



n 2(l + 6AVc)w, 



^fc - — ^2 (Bl) 



Here Sk is the function Sat (A), with N replaced by k and A is fixed. Let us stress the relation between the N- 
dimensional integral and the function Fj^: 



1 



The aim of the Gelfand-Yaglom construction is to find in the continuum limit such a differential equation in variable 
T ~ A;. A that its solution is connected to the continuum functional integral by relation: 
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In the spirit of the Gelfand-Yaglom construction we are going to express the relation for Fk+i by help of Fk and 
We have: 



n 2(1 + 6AVc)cj, 2(1 + 6AVcK+i J] 2(1 + bA^/c)u 



Fi 



fc+i 



i=0 



1=0 



(B2) 



2(1 + 6A7c) n 2(1 + bl\^/c)uji n 2(1 + 6A7c) 



Regarding to the definition of the function Fk, in Eq. (jBip . we have: 



Ffe+i = 2(1 + bAyc)Fk ^ - Ffe_i 



02 



After some algebra we find: 



Fk+i — 2Fk + Fk-i — 2{Fk — Fk-i) 



C2 _ c2 
c2 

•^k+l 



2&A2 



-Fk — Fk-i 



'^k+1 



(B3) 



-Fk 



^k+i 



We need not to know the structure of functions Sk to derive the identities: 



^k ^ ^k+l — ^"^Sk+liSk+l ~ Sk) + {Sk+l ~ Sk) 



and 



Sl_i — 2Sl + Sl_^_l — 2{Sk+i — Sk)"^ + ^Sk+iiSk+i — 2Sk + Sk^i) — 

~ A{Sk+i - Sk){Sk+i - 2Sk + Sk-i) + {Sk+i - 2Sk + Sk-i)^ . 
Inserting these identities into Eq. (jB4p we find a difference equation which, divided by A^, takes the form: 

oL / c^Sk+1 — 2Sk + Sk-i „ / Sk+i — Sk 
2o/c — 2 . „ 2 



(B4) 

(B5) 
(B6) 



Fk+i — 2Fk + Fk-i , .Fk — Fk-i Sk+i — Sk „ 
- 4 : : — = I'k 



A2 
2/ 



A ASk 



A-'Sk 



AOi + A'Oi , 



ASk 



(B7) 



where 



4b fSk+i_Sk\^Jjk^Fk-i] 



c V '^Sk+i 



4Fk-i 



f Sk+i — Sk\ ( S'fc+i — 25*^ + S'fc-i 



V ASk + 1 y V A25fe+i 



1-^2 = — i'k — -r-^ — fk-\ 



ASk+i 



A^Sk, 



In the contimmm limit A ^ we replace k.A by r. Under the condition 

hm^(A d + A^ O2) = , 



(B8) 



(B9) 



we obtain the differential equation: 



9^Fir)+i-^Fir)^lnSir)^Fir) 



^-2|^ln.(.)-4(^ln.(.) 



(BIO) 
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The initial conditions for S{t) continuous and finite in t = are: 

m ^ 



^2(0)' 

— F(0 = lim 



(Bll) 



52(0) 



APPENDIX C: EVALUATION OF THE LEADING PART OF Sn 

In this appendix we evaluate the finite range ki summations by the recurrence method. Let us start with summation 
over the index fcjv-i of the Eq.([5]). The finite sum to be done is: 



Ko 



zr = E 



1=0 



1 



v/2^(l + 6A2/c) (2fcAr_i)! 



r(A:jv-2 + fcw-i + l/2)X'_fc„„2_fc„_i_i/2 (z) 



1 



r(fcjv-i + l/2)I?_fc„_,_i/2 (zo) 



V27r(l + feA2/c)cJo 

Let us recall the dependence of the function V on the parabolic cylinder function D: 

and also the definitions of the variables ^, ujq, z and zq: 

c(l + bA2/c) _ c(1/2 + &A2/c) 



(CI) 



v/2aA3 
1 



LUq 



v/2aA3 
1/2 + 6A2/C 



l + fcA2/c' 1 + 6A2/C ■ 

We use the asymptotic Poincare-type expansion of the parabolic cylinder function, which for V means: 

I?_,.„_V2(^o) ^ zo'"^'^' e^'/' ^-fe„-i/2(^o) = E(-l)^"^^^^^^-^+£j(fc„,zo) 



(C2) 



In the last relation, denotes the number of terms of the asymptotic expansions convenient to take into account. 
We apply asymptotic expansion for the function I?_fc„_j_i/2 (zo) in Eq. ljCip . In the truncated sum we interchange 
the order of the finite summations over indices kw-i and j. We replace V by D, therefore a corresponding power of 
the variable z will play an important role. In this way we obtain the relation: 



^cut 



r(i/2) 



E 



{-ly exp(z2/4) 



^fc„_.2 + l/2 



v/27r(l + feA2/c)t^o i! (24)^' ^/M^ + ^A/c) 



(C3) 



E 



(fcw-i)! 



r(A:jv-i + kN-2 + 1/2) (few-i + 1/2)2^. i:'-fc„_i-fc„_,_i/2 (z) 



fcjv-i=0 

where we simplified the calculations by identities: 

(2fc)! = 2^'= A:! (1/2)^, r(fc + 1/2) = r(l/2)(l/2)^. 

Let us study in detail the sum 



i Uiv-i 

{kN-l)\ 



r(/cAr_i + fc7V-2 + 1/2) (few-i + 1/2)2^. D-fc„_i_fc„_,-l/2 (z) 



(C4) 
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This sum is uniformly convergent, therefore we can extend the summation up to infinity by adding the corresponding 
terms, which appear also in the remainder with opposite sign. To be able to provide the sum over the index kj^^i, 
we must modify the Pochhammer symbol 

(fcjv-i + 1/2)2^. = (fcjv-i + 1/2). • • • .(fcAT-i + 1/2 + 2j - 1) . 

We see, that this object is a polynomial in the variable kjy^i of the 2j — th order. We rewrite the polynomial in 
another form: 

min (2j,feiv-i) 

(fcAT-i + 1/2)2^. = 2^ a/ 



The coefficients a^-' are given by recurrence procedure from the relation: 

s ' XT' ^ s s 

kN-i=Q ^ ' 2=0 fcjv-i=j ^ " ^ ' 

From the above definition, we find the recurrence equation: 

o}l = {k-\|2 + ^)a\-^ ^a\zl (C6) 

while the initial conditions are: 

= 1 , = (1/2)^ , a^+i = 



The solution of this recurrence equation is: 

(1/2) 



(C7) 



* ; (1/2), 

Inserting all these replacements into Eq. (|C4p . with help of the identity 

r(fcAr-2 + + 1/2) = r(fcjv-2 + i + 1/2) (A:w„2 + « + l/2)fc„_^_, 

after some algebra, introducing a new summation index k = kj^^i — z, we obtain the formula: 

X rikN-2 + I + 1/2) j ^^^g- (fcjV-2 + * + l/2)fc D-fe„_,-._l/2-fc (2) 

In the above relation we extend summation over the index k up to infinity. The sum over k is now prepared for 
application of the identity: 



fc=0 

The result of the first recurrence step, replacing D by P, reads: 



g.V4 ^ D^,^,{x) = e(--*)^/4 Z?_,. (:r - t) (C9) 



^r*- , ^ , ^"^'^ y T"!' ( ^^V^(fcw_2 + ^ + l/2)2?_. ,_,_i/2 (^i) 

(CIO) 

where 

^1 = ^ (1 - CV(4^o)) , c^i = ^ = 1 - eV(4wo) . 



17 



For the following summation over the index A;jv-2 we have: 

1 1 



V2(l + 6A2/c)a;o ^2(1 + bA^/c)uji 



(Cll) 



X E 

feN-2=0 



j=0 -I- ^'^^ I i=o 



V27r(l + 6A2/c) (2A;jv-2)! 

r E ( ^ 



■r(fcAr_3 + A;jv-2 + l/2)r'_fe 



JV-3 — ^iV-2- 



1/2 (2;) 



4a;oWi 



r(fcjv-2 + i + l/2)2?_fc^_,_i_i/2 (^1) 



where 



(C12) 



will define the first step of the new recurrence relation. 

Now, due to the uniform convergence of the sum over the index kM-2 we will evaluate the leading part of Z2 as 
a finite sum over index kM-2- In the finite sum, we use the asymptotic Poincare-type expansion of the parabolic 
cylinder function 'D_^.^_^_^_l/2 {zij. We have then in the relation for ^2"* finite summations only and we change 
the order of the sums. We have: 



crcut 

^2 



1 



1 



1 



^2(l + 5A2/c)a;o ^2(1 + 6A2/c)u;i yj2TT(l + bA^/c) 

,2 \i J' 



h^'-^'^'^h^ ' V-oc.rj — 



X E i^TTI r(fcjv_3 + fcjV-2 + 1/2) £'-fejv-3-feiV-2-l/2 (2;) (fcjV-2+l/2)2i+i 

feiV-2=0 

where A = ^/2. Summing over the index kN-2 as in the first recurrence step, we have: 



(C13) 



■^2 



1 



1 



(^2) 



EE ^ 



i=0 



where one defines new variables: 



V2(l + 6A2/c)a;o v'2(l + 6A2/c)a;i ^2(1 + 6A2/c)a;2 ^ j!;!(2z2)j+i ^2; 

2j / /12 \ ' 2;+s 

E 

p=0 



-21 (C14) 



■' / /12 \' / 42 \P 

EWfi^j E«rM^) r(fc;v-3 + P + i/2p. 



-feiV-3-p-l 



72(^2) , 



Z2 = ZUJ2, 0^2 = 1 



The summations over indices j, I is done as follows: 



J J 



j=0 1=0 

J 



(C15) 



M=0 



3=0 



2J 



J 



IJ.=J+1 



(2^2)/^ 



9{j) HP' - j) 
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The first term on the right hand side of the above relation will contribute to the leading part of Eq. (jC14[) , while the 
second term, where index [i> J ^ will contribute to the remainder, due to the term z~'^^ ~ which may be made 

as small as possible by choosing J" properly. Interchanging the order of summations: 

2j 2fj.-2j+i IJ. 2j 

E E E 

i=0 p=a p=a i=max [0, p-2/i+2j] 

we find the result of the second recurrence step: 

' ^2(1 + 6A2/c)ljo y/2{l + 6A2/c)cji ^2(1 + bAyc)uj2 

where the second recurrence step of the function (2)^^ is defined by: 

/ \ 1 2j / a2 



2j 2f,-2j+i 



We can see that after A recurrence steps the result of the A summations over the indices ki can be read@: 



(C17) 



^T* = Ij] , ^ =1 (^a)-''"-^-^ (C18) 

X E «(^)^(fc.-A-r+P+l/2)P-._-,-V2(.A). 

We have evaluated the recurrence relations: 

l/2 + 6AVc , e 
WA-i l + oA-^/c 2 

and introduced the recurrence definition for the function (A)^'^: 



j=0 ^ ^ ^A-l i=max[0,p-2/.+2j] \'^A-2WA-l/ 



(C19) 



where the recurrence procedure begins from (1)^"' given in Eq. (jC12[) . 

After the last recurrence step, for A = iV — 1, we are left with the relation of the form (jC18P where fcjv_A-i = 
and the index of the V function is only —p — 1/2. We expand the D^p^i/2{zN-i) as in all previous recurrence steps 
and find the relation: 

Zn-1 = \ n -7^7TTtxTjt= r E ^ ^ 



l\ V2(l + 6AVc)w. J t'o 



, ^ / .. \ 1 2i / ^2 

X 

Following the definition of the aj symbols, we have: 

(1/2).,+. = 



19 



Then in the last part of the preceding equation we read: 

E f ) ^ E ( - i)f «r - (C21) 

Following the calculations done in this Appendix, we conclude, that it is possible to realize summations in the exact 
formula for the N— dimensional integral at least by help of asymptotic expansions of the parabolic cylinder functions. 
It is possible to provide the continuum limit of our result and there are no additional terms contributing to the result 
in the continuum limit. The result reads 

This expression is sufficient for the calculation of the continuum unconditional Wiener measure functional integral 
by the Gelfand-Yaglom procedure leading to the differential equation of the second order. The second part of the 
relation (jC22[) represents the expansion of an unknown function. 

In the following calculation the key role play objects LOi defined by the recurrence relation 

U)i = 1 - 



with the first term 

Wo = 1/2 + B, 

where 

6AVc 



B = 



2(l + 6A2/c)' 



Ui defined in this way are represented by continued fractions. The continued fraction can be represented by a simpler 
relation as the solution of its n — th convergent problem of continued fractions [l5l |. Let us shortly explain this 
procedure. 

Let us have a continued fraction of the form: 

u! — ai + - 



«2 + ^ 



The n — th convergent is defined as 



where p„ and q„ are given by equations: 



Pn 
= — 



Pn = an Pn-1 + K Pn~2 
In = In Qn-l + b„ qn-2 

(C23) 



Solutions of these recurrence equations have the form: 



Pn = WiPi + W2P2 
Qn = WlPi + W2P2 

(C24) 



where pi^2 are solutions of the characteristic equation, in our case a homogenous one: 

- ttnP - 6„ = 



For the continued fraction in question we have: 



20 



and the solution of the characteristic equation is: 



-(1± -4/12 ) 



Pl,2 2 

The constants Wi and Wi are fixed by lvq and uji terms, that adjust the initial conditions: 

Po = 1 + 2B 

Pi = l + 2B-A^ 

90-2 

qi = 1 + 2B 

The n — th convergent method solution is completed by the relations: 

1 r.. / / 2B 

^1.2 = 2 



(1 + 25) ±(71311^4 



2B 

Wi2 = 1 ± 



VI - 4yl2 ' 

the particular characteristic follows from the above solution: 

Pn = gn+i 

which simplifies our calculation significantly. In forthcoming calculations we have introduced more convenient vari- 
ables: 

and also performed the replacement: 



^2 n2 



In what follows, we: 

- replace the summation index i by the summation index j defined by i = 2A — j . 

- interchange the order of summations over indexes j and A. 

- rewrite the recurrence relation Eq. (|C19[) as follows: 

= il iAnn~\v-, E (^Qa-2Qa-i)''-^' (A - 1)^^^- { ^ \ {Ql^.T-' ,pe<0,2^,>, (C25) 

the recurrence starts from the term (jC12|) 



APPENDIX D: EVALUATION OF THE REMAINDER TO THE LEADING PART OF 5jv 

Let us describe the evaluation of the remainder to the leading part of the sum calculated in Appendix C. The 
remainder consists of the infinite sum of the original series, the finite sum over the remainder of the Poincare expansion 
of the parabolic cylinder function, the second term of the finite sum Eq. (|C15p and the relation added to the leading 
part for possibility to perform an infinite summation over the parabolic cylinder function: 

R{J, Ko)^~J2 -uo 2)j E tItv ^(^-1 + + ^/mh + h+1 + 2j + l/2)P_,.,_,_fe,_i/2 (z) + 

j=0 ^ ^ ) k, = Ka ^ 

+ E 7liTT^(^'-i + + l/2)r(fc, + k,+, + l/2)I?_fc,_^_fc,_i/2 (z) ej{k, + /c,+i, z) + (Dl) 

ki=0 ^ 

+ E Tfiv^^^*"' ^ ^ + ^^+1 + l/2p-fe._i-fe.-l/2 (^)2?-fe.-fc.+ i-l/2 (z)- 

k,=Ko + l ^ 
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The finite sum in remainder: 



J2 tIf^^^^'-i + + V2)r(fc. + h+i + l/2)I?_fe,_,_fe._i/2 (z) ej{h + h+i,z) 

ki=Q ^ 



is bounded by the relation (|26p : 
M 



tItv ^(^'-1 + + V2)r(fc. + h+i +2J + l/2)P_fc,_^_fc^_i/2 (z) 



where 



{J-l)\ {2z^y {2h) 



2z^ ( 4(5 a? 

7W = max ( -i i^2(^/2, 1/2; J/2 + 1; 1 - — ) exp iF2{ll2, 1/2; 3/2; 1 - ^) 



z2 - 2a 



22 - 2a 



(D2) 



(D3) 



(D4) 



and a = fci + /j^+i, = 1,2,..., iiTo • Let us remember that ^ . Since we have the freedom to choose the 
parameter J, the upper bound on this contribution to the remainder can be made as small as necessary power of N . 

In the asymptotic region of ki > Kq we expand one of the function V to double asymptotic expansions proposed 
by Temme [16| : 



.fe=o 



where the following quantities were introduced: 

1 



A = — , Wo = 77 [Vl + 4A — 1] , ^ = - + — A — A In Wo + A In A , 
2 z 

the functions /fc(A) are calculated in [16]. We find for the infinite part of the sum decompositions: 



(D5) 



(0 



,5.r (2^^)' 



r(/c,_i + h + l/2)T{h + h+i + l/2)2?_fe,_^_fe,_i/2 (z)I?-fe.-fc,^,„i/2 (^) 



E 



i^r(fc,_i + fc, + i/2)r(fc, + fc,+i + i/2)2?_,._,_,,_i/2 (^) (7+ 4^^1/1 

exp (— ylz^) 



E 

.fc=0 



A (A) 



+ 



(0 



ki=K, 



- ^kcf^^'-^ ^ + i/2)r(fc, + fc,+i + i/2)i)_fe._,_fe._i/2 (^)-7fqr^^ 

+ 1 



— i?„(a, z) 



(D6) 



It was shown in [Sj] that the last part of this contribution can be made as small as we need due to the freedom in the 
choice of the parameter n representing the number of the functions /fe(A) taken into account in the Temme double 
asymptotic decomposition ^\&\ of the parabolic cylinder function. 

Now we are going to estimate the last part of the remainder, corresponding to the difference of the series: 



ki=Ka + 
J 



E 

.fc=0 



/fc(A) 



E 



E 



(0 



Iki 



r(fc,_i + fc, + i/2)r(fc, + fc,+i + 2j + i/2)i?_fe,_^_fc^_i/2 (z) 



3=0 ■' ^ ' k,=Ko + l ^ ' 

Since both series are uniformly convergent, we exchange the order of summations, to obtain: 



(D7) 



2fei 



E iliv^^^'-' + + l/2)r(fc, + A:,;+i + l/2)P_fc,_^_fe^_i/2 (^) 



exp {—Az^) 
(l + 4A)i/4 



E 

,fc=0 



/fe(A) 



,2fe 



^ (-lP(fc. + fc.+l + 1/2)2, 



j! (2z2y 



(D8) 
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The double asymptotic expansion reduces to Poincare - type expansion [8| when a is fixed after expanding the 
quantities in A = a/ for small values of this parameter. In the difference 



exp {—Az'^) 
(l + 4A)i/4 



LA.— 



f (-lV(fe. + fcm + l/2),, . 



all terms where fc,j < min(n,i/) cancel one another and the rest of terms is proportional or smaller a 

As in the case of previous contributions to the remainder, this means that this part of the remainder can be made as 

small as we need in the power of 



APPENDIX E: DEFINITION OF THE MATRICES 



The definition of the matrices A''(A - 1), C''(A - 1), M''(A - 1) is the following: 

1. The A'' is the lower triangular matrix with zeros over the main diagonal of the dimension (2/i + l)(2/x + 1). The 
principal minor of the dimension (2d + \){2d + 1) is non-zero only with elements: 



{A'^(A-l)} 



a. 



2d~j 
2d-p 



2. The C^_;^(A) is the upper triangular matrix with zeros under the main diagonal of the dimension (2/i+ + 1). 
The nonzero elements form the main minor of the dimension {2d + l){d + 1) with A*'' column: 

where p = 0, 1, 2A and A = 0, 1, d . 

3. The f. is the upper triangular matrix with zeros under the main diagonal of the dimension (/i + + 1). 
The nonzero elements form the main minor of the dimension (d + l){d + I) : 



Ml,(A - 1) = ( A ) (Qi-i)"""' d>k>X>0. 
For a product of two lower-triangular matrices A'^{2) and A^^(l) we haveg: 



P P „2/3-i ^2^2 -A 

i:{A'-(2,}„{A'.,i)},,,^ 1:^35^(3^^ (El, 

. ,-2f.-A) (4i.-2A)! [^..-2..^ /^P-AW V'' ( ^ 

{4h-2p)l{p-Xy. ' y ' )\AQ^Q2) \AQ2Q 

where e is an auxiliary variable. 

We have used two identities for the summation over the index j and the definition for index Ij 



and 



2/3-. _ ,-2(p-A) (4/2 -2A)! / P - A \ (4/3 -2j-)! 

«2/3-P «2/.-, (4^^ -2p)\(p-\)\\3-\) (4/2 - 23)\ 



(4/3 2j)! _ p,il3-U2(Al3~2o\\ 

(4/2-2J)! -"^^ ^1^=^ 



Ij = d - {ij H hu) 
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In the above relation the summation over the index j can be performed exphcitly. Introducing a new auxiHary 
variable: 



we find: 



2(p- 



-A) (4/2 ^„ 



(4/3-2p)!(p-A)! « ]^p-2h \AQ^Q2 AQ^Q 



p-X' 



where — 3/49| + 3^9| + is calculated from 9^ . 
For the resulting product of all matrices A^(fe) we find: 



{A''-''^(A - 1) * A''-''^-*^-i(A - 2) * • • • * A'^-^'^-^"-^ ''(l)}p A = 



^ g-2(p-A) (4/2 -2A)! 4(^^+i3 + ...+^^)^ 

(4/A-2p)!(p-A)! 



,m=2 



6 



62 ■ ■ -Ca-i 



^p-2/„+i ^^Q^Ql ^Q3Q2 AQaQ 



A-1 



p-A 



(ai/Cm — 1) 



By symbol * we indicate the product of matrices. Evaluating the product of matrices: 



jiyjd-^A-iA-i i2 (1) ^ . . . ^ Md-»A-iA-i (A - 2) * M'^-*'^ (A - 1)| 

we use that M^^+i(j) is a one-column matrix with Ij+i non-zero elements in the (j -|- 1) column: 



A 



where A = 0, 1, • • • , Ij+i- Product of such matrices is a one-column matrix with elements: 



vd—iA—iA- 



>-{l). 



A, /a 



A j U2 



/a a r)4(^2-A)^4(/3-/2) ^4(/a-/a-i) 
J^_^ 1^1 ^2 '''^A-1 



From the definition pop of recurrence steps we have for the matrix C^^(l) nonzero elements: 

{C^^l)}, 



4A 

2A 

"2A-J 



with conditions for indexes: 



< j < 2A < 2/2 < 2/i 



Collecting all partial results together, inserting them into Eq. (|3ip and remembering that for function 5a defined 
in Eq. ([29|) only matrix elements {C(A)^^}2^ are important, we find the result (|32|) . 
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